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ABSTRACT 

The non-Gaussian cold spot detected in wavelet space in the WMAP 1-year data, is detected again 
in the coadded WMAP 3-year data at the same position (b — —57°, I — 209°) and size in the sky 
(« 10°). The present analysis is based on several statistical methods: kurtosis, maximum absolute 
temperature, number of pixels below a given threshold, volume and Higher Criticism. All these 
methods detect deviations from Gaussianity in the 3-year data set at a slightly higher confidence 
level than in the WMAP 1-year data. These small differences are mainly due to the new foreground 
reduction technique and not to the reduction of the noise level, which is negligible at the scale of 
the spot. In order to avoid a posteriori analyses, we recalculate for the WMAP 3-year data the 
significance of the deviation in the kurtosis. The skewness and kurtosis tests were the first tests 
performed with wavelets for the WMAP data. We obtain that the probability of finding an at least 
as high deviation in Gaussian simulations is 1.85%. The frequency dependence of the spot is shown 
to be extremely flat. Galactic foreground emissions are not likely to be responsible for the detected 
deviation from Gaussianity. 

Subject headings: methods: data analysis - cosmic microwave background 



1. INTRODUCTION 

The Cosmic Microwave Background (CMB) is at the 
moment the most useful tool in the study of the ori- 
gin of the universe. A precise knowledge of its power 
spectrum constrains significantly the values of the cos- 
mological parameters which determine the cosmological 
model. The 1-year Wilkinson Microwave Anisotropy 
Probe data (WMAP, Bennett et al. 2003a), measured 
the anisotropies of the CMB with unprecedented accu- 
racy, finding that the standard model fits these data. 
A flat A-dominated Cold Dark Matter (ACDM) uni- 
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verse with standard inflation explains most of the ob- 
servations confirming the widely accepted concordance 
model. According to standard inflation, the tempera- 
ture anisotropies of the CMB are predicted to represent 
a homogeneous and isotropic Gaussian random field on 
the sky. A first Gaussianity analysis found the data to 
be compatible with Gaussianity (Komatsu et al. 2003). 

Several non-Gaussian signatures or asymmetries were 
detected in the 1-year WMAP data in subsequent works. 
A variety of methods were used and applied in real, har- 
monic and wavelet space: low multipolc alignment statis- 
tics (de Oliveira-Costa et al. 2004, Copi et al. 2004, 
2005, Schwarz et al. 2004, Land & Magueijo 2005a,b,c, 
Bielewicz et al. 2005, Slosar & Seljak 2004); phase cor- 
relations (Chiang et al. 2003, Coles et al. 2004); hot 
and cold spot analysis (Larson & Wandelt 2004, 2005); 
local curvature methods (Hansen et al. 2004, Cabella et 
al. 2005); correlation functions (Eriksen et al. 2004a, 
2005, Tojeiro et al. 2005); structure alignment statistics 
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(Wiaux et al. 2006); multivariate analysis (Dineen & 
Coles 2005); Minkowski functionals (Park 2004, Eriksen 
et al. 2004b); gradient and dispersion analyses (Chyzy et 
al. 2005); and several statistics applied in wavelet space 
(Vielva et al. 2004, Mukherjee & Wang 2004, Cruz et 
al. 2005, 2006, McEwen et al. 2005a and Cayon, Jin & 
Treaster 2005). 

The recently released 3-year WMAP data with higher 
signal to noise ratio is key to confirm or disprove all these 
results. 

In the 3-year papers, the WMAP team (Hinshaw et al. 
2006) re-evaluates potential sources of systematic errors 
and concludes that the 3-year maps are consistent with 
the 1-year maps. The exhaustive polarization analysis 
enhances the confidence on the accuracy of the temper- 
ature maps. The ACDM model continues to provide the 
best fit to the data. 

Spergel et al. (2006) perform a Gaussianity analysis 
of the 3-year data. No departure from Gaussianity is 
detected based on the one point distribution function, 
Minkowski functionals, the bispectrum and the trispec- 
trum of the maps. The authors do not re-evaluate the 
other statistics showing asymmetries or non-Gaussian 
signatures in the 1-year data. 

The aim of this paper is to check the results of Vielva 
et al. (2004), Cruz et al. (2005), Cayon, Jin & Treaster 
(2005) and Cruz et al. (2006), (hereafter V04, C05, CJT 
and C06 respectively) with the recently released WMAP 
data. All these analyses were based on wavelet space. 
In particular the data were convolved with the Spherical 
Mexican Hat Wavelet (SMHW). Convolution of a CMB 
map with the SMHW at a particular wavelet scale in- 
creases the signal to noise ratio at that scale. Moreover, 
the spatial location of the different features of a map is 
preserved. 

V04 detected an excess of kurtosis in the 1-year 
WMAP data compared to 10000 Gaussian simulations. 
This excess occurred at wavelet scales around 5° (an- 
gular size in the sky of m 10°). The excess was found 
to be localized in the southern Galactic hemisphere. A 
very cold spot, called the Spot, at galactic coordinates 
(b = — 57°, I = 209°), was pointed out as the possible 
source of this deviation. 

C05 showed that indeed the Spot was responsible for 
the detection. The number of cold pixels below sev- 
eral thresholds (cold Area) of the Spot was unusually 
high compared to the spots appearing in the simulations. 
Compatibility with Gaussianity was found when mask- 
ing this spot in the data. The minimum temperature of 
the Spot was as well highly significant. 

C06 confirmed the robustness of the detection and 
analysed the morphology and the foreground contribu- 
tion to the Spot. The Spot appeared statistically robust 
in all the performed tests, being the probability of find- 
ing a similar or bigger spot in the Gaussian simulations 
less than 1%. The shape of the Spot was shown to be 
roughly circular, using Elliptical Mexican Hat Wavelets 
on the sphere. Moreover the foreground contribution in 
the region of the Spot was found to be very low. The 
Spot remained highly significant independently of the 
used foreground reduction technique. In addition the 
frequency dependence of the Spot was shown to be ex- 
tremely flat. Even considering large errors in the fore- 
ground estimation it was not possible to explain the non- 



Gaussian properties of the Spot. 

CJT applied Higher Criticism statistics (hereafter HC) 
to the 1-year maps after convolving them with the 
SMHW. This method provided a direct detection of the 
Spot. The HC values appeared to be higher than 99% of 
the Gaussian simulations. 

Note that although the Spot has not been detected 
in real space, this structure exists but is hidden by 
structures at different scales. The convolution with the 
SMHW at the appropriate scale, amplifies the Spot, mak- 
ing it more prominent. 

Several attempts have been made in order to explain 
the non-Gaussian nature of this cold spot. Tomita (2005) 
suggested that local second-order gravitational effects 
could produce the Spot. Inoue & Silk (2006) consid- 
ered the possibility of explaining the Spot and other large 
scale anomalies by local compensated voids. Jaffe et al. 
(2005a) and Cayon et al. (2006) assumed an anisotropic 
Bianchi VII/j model showing that it could explain the ex- 
cess of kurtosis and the HC detection as well as several 
large scale anomalies. On the other hand, McEwen et al. 
(2005b) still detect non-Gaussianity in the Bianchi cor- 
rected maps. Jaffe et al. (2005b) proved the incompati- 
bility of the extended Bianchi models including the dark 
energy term with the 1-year data. Adler et al. (2006), 
developed a finite cosmology model which would explain 
the Spot and the low multipolcs in the angular power 
spectrum. Up to date there are no further evidences of 
the validity of any of the above suggested explanations. 

Our paper is organized as follows. We discuss the 
changes in the new WMAP data release and the simula- 
tions in section §2. The analysis using all the mentioned 
estimators is described in section §3. In section §4, the 
significance of our findings is discussed. We analyse the 
frequency dependence of the Spot in section §5, and our 
discussion and conclusions are presented in sections §6 
and §7. 

2. WMAP 3- YEAR DATA AND SIMULATIONS 

The WMAP data are provided at five frequency-bands, 
namely K-band (22.8 GHz, one receiver), Ka-band (33.0 
GHz, one receiver), Q-band (40.7 GHz, two receivers), 
V-band (60.8 GHz, two receivers) and W-band (93.5 
GHz, four receivers). Foreground cleaned maps for 
the Q, V, and W channels are also available at the 
Legacy Archive for Microwave BAckground Data Anal- 
ysis (LAMBDA) web site 3 . 

Most of the 1-year Gaussianity analyses were per- 
formed using the WMAP combined, foreground cleaned 
Q-V-W map (hereafter WCM; see Bennett et al. 2003a). 
CMB is the dominant signal at these bands and noise 
properties are well defined for this map. The de-biased 
Internal Linear Combination map, (DILC) proposed by 
the WMAP team, estimates the CMB on the whole sky. 
However its noise properties are complicated and regions 
close to the Galactic plane will be highly contaminated 
by foregrounds. Chiang, Naselsky & Coles (2006) find 
evidences for the foreground contamination of the DILC. 
Therefore we will still use the more reliable WCM in the 
3-year data analysis. 

Hinshaw et al. (2006) describe some changes in the 
3-year temperature analysis with respect to the 1-year 

3 http://lambda.gsfc.nasa.gov 
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Fig. 1. — Image showing an azimuthal projection of a 22° X 22° patch from the WCM HEALPix map with resolution nside = 256, 
centered on the Spot and in fiK. In the first row we have the 1— year and 3— year images of the Spot in real space, whereas in the second row 
the Spot is shown at wavelet scale Rg. The image is divided in 1024 X 1024 pixels and the y-axis is oriented in the Galactic north-south 
direction. 



one. Coadding the three years of observations reduces 
the instrumental noise. The 3-year maps have w 3 times 
lower variance. Refinements in gain calibration and beam 
response models have been implemented and a new fore- 
ground reduction technique has been used. The latter 
seems to provide a better correction than the one ap- 
plied to the first year data. As discussed in C06 the 
Galactic foreground estimation is a very important issue 
in Gaussianity analyses. The exclusion masks defined by 
Bennett et al. (2003b) have not been modified, except 
for the inclusion of 81 new point sources in the kpO mask. 
This mask excludes the highly contaminated pixels close 
to the Galactic plane. 

Despite these changes the 3-year maps have been found 
to be consistent with the 1-year maps by the WMAP 
team. 



V04 and C05 performed a very careful analysis in or- 
der to study the power spectrum and noise dependence 
of the kurtosis and cold Area estimators. Considering 
different power spectra within the la error band of the 
1-year data, the differences in the significance of the kur- 
tosis were found to be negligible (see Figure 11 in V04). 
The Area of a particular spot was neither affected by 
the power spectrum (see section 5.3 in C05). The results 
were almost noise independent. The convolution with 
the SMHW reduces considerably the noise contribution. 
Even if similar results are expected, we perform 10000 
Gaussian simulations of the 3-year coadded data follow- 
ing the same steps as for the 1-year simulations. The 
only differences between the 3-year and the 1-year sim- 
ulations are a lower noise contribution and a very slight 
variation in the power spectrum used to generate the 
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Fig. 2. — WCM kurtosis values for the 1— year (asterisks) and the 
3-year data (circles). The acceptance intervals for the 32% (inner), 
5% (middle) and 1% (outer) significance levels, given by the 10000 
simulations are also plotted. 



simulations. For a detailed description of the simulation 
pipeline, see section 2 of V04. 

We will use all these maps in the HEALPix pixclisation 
scheme (Gorski et al. 2005) 4 with resolution parameter 
Nside = 256. 

3. ANALYSIS 

Our aim in this section is to repeat the same tests per- 
formed in V04, C05, CJT and C06 but with the 3-year 
data. Then we will compare the new results to the old 
ones. One can see the region of the Spot in real and 
wavelet space at scale 5° for both releases of the WMAP 
data in Figure^ In real space the 3-year data image ap- 
pears clearly less noisy, whereas the wavelet space images 
present only very small differences. 

In V04, data and simulations were convolved with the 
SMHW at 15 scales, namely (R t = 13.7, R 2 = 25, 
i? 3 = 50, i? 4 = 75, R 5 = 100, i? 6 = 150, R 7 = 200, 
R s = 250, R g = 300, i?m = 400, R n = 500, R u = 600, 
i?i 3 = 750, i?M = 900 and i?i 5 = 1050 arcmin). The 
SMHW optimally enhances some non-Gaussian signa- 
tures on the sphere (Martinez-Gonzalez ct al. 2002) and 
has the following expression: 
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is a normalisation constant: N(R) = 
R 2 /2 + i? 4 /4. The distance y on the tangent 
plane is related to the polar angle (9) as: y = 2tanf9/2. 

We will use the same 15 scales in our present analy- 
sis, considering those estimators where non-Gaussianity 
was found in the 1-year data, namely kurtosis, Area, 

4 http:/ /www. eso.org/science/healpix/ 



Max, HC and a new one, the volume. The definitions of 
each estimator will be given in the following subsections. 
Analyses were also performed in real space, which will 
be referred as wavelet scale zero. In real space, the data 
are found to be compatible with Gaussian predictions 

In the following subsections we will give the upper tail 
probabilities of the data at one particular scale. The up- 
per tail probability is the probability that the relevant 
statistic takes a value at least as large as the one ob- 
served, when the null hypothesis is true. 

In section §4 we will give a more rigorous measure of 
the significance, considering the total number of per- 
formed tests to calculate the p-value of the Spot. The 
p- value is the probability that the relevant statistic takes 
a value at least as extreme as the one observed, when the 
null hypothesis is true. In our case, the null hypothesis 
is the Gaussianity of the temperature fluctuations. 

3.1. Kurtosis 

Given a random variable X, the kurtosis k is defined 
as k{X) — p^pfgTp — 3. In V04 the kurtosis of the 
wavelet coefficients was compared to the acceptance in- 
tervals given by the simulations. In FigureElthe kurtosis 
of the 1-year data are represented by asterisks and the 3- 
year data by circles. Hereafter we will use these symbols 
to represent 1-year and 3-year data. Both are plotted 
versus the 15 wavelet scales. Scale corresponds to real 
space. The acceptance intervals given by the simulations 
will be plotted in the same way in all figures: the 32% 
interval corresponds to the inner band, the 5% interval to 
the middle band and the 1% acceptance interval, to the 
outer one. As expected, the acceptance intervals remain 
almost unchanged with respect to those obtained from 
1-year simulations. This will happen as well for all the 
other estimators. The 3-year kurtosis values follow the 
same pattern as the 1-year ones, confirming the initial 
results. However there are slight differences at the scales 
where the deviation is detected, being the kurtosis even 
higher in the 3-year data. The most significant deviation 
from the Gaussian values, occurs at scale Rg = 5°. In Ta- 
ble ^ we list the kurtosis values at scale Rg, considering 
the 1-year data as published in 2003, the 1-year data re- 
lease applying the changes in the data analysis described 
in Hinshaw et al.(2006), and the coadded 3-year data. 
The biggest difference is found between both releases of 
the 1-year data. The kurtosis value of the 1-year data 
increases « 7%. This may be due to the new foreground 
reduction technique. As expected the noise reduction 
due to coadding the three years of observations, implies 
a much lower increase in the kurtosis, since the noise con- 
tribution in wavelet space is very small. The upper tail 
probabilities (i.e. the probabilities of obtaining higher 
or equal values assuming the Gaussian hypothesis) are 
given in the right column of Table ^ Hereafter we will 
compare the first release of the 1-year data with the 3- 
year data. 

Analysing both Galactic hemispheres separately, we 
obtain the results presented in Figure El Again the kur- 
tosis follows the same pattern as in the 1-year results. 
As expected, the deviation appears only in the southern 
hemisphere and it is slightly higher in the 3-year data. 
The upper tail probability obtained in V04 was 0.11% 
at scale Rj in the southern hemisphere, whereas now we 
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Fig. 3. — As in Figure 121 but for the northern (left plot) and 



TABLE 1 

KURTOSIS VALUES AT SCALE R$ 



Data 


kurtosis 


probability 3 


1-year data (2003) 


0.836 


0.38% 


1-year data (2006) 


0.895 


0.28% 


3— year data 


0.915 


0.23% 



Note. — Kurtosis values of different WCM versions at scale 
Rg. The right column gives the probability of obtaining a higher 
or equal value in Gaussian simulations. 



have 0.08% again at scale Rj. The deviation from Gaus- 
sianity is localised in the southern hemisphere because 
the Spot is responsible for it (see C05). 

3.2. Maximum statistic 

Given n individual observations Xi , Max is defined as 
the largest (absolute) observation : 

Max n = max{\Xi\, \X 2 \, \X n \}. 

The very cold minimum temperature of the Spot, was 
shown to deviate from the Gaussian behaviour in V04. 
In this work and in C05, C06 the minimum tempera- 
ture estimator was used to characterise the Spot whereas 
in CJT the chosen estimator was Max. As Max is a 
classical and more conservative estimator, we will use it 
in the present paper instead of the minimum tempera- 
ture. Our n observations correspond to values in real or 
wavelet space (normalized to zero mean and dispersion 
one). The Spot appears to be the maximum absolute 
observation of the data at scales between 200 and 400 
arcmin. In Figure^] the 1-year and 3-year WMAP data 
values of Max are compared to those obtained from the 
simulations. As for the kurtosis, both data releases show 
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Fig. 4. — Maximum absolute observation versus the 15 wavelet 
scales. Again the circles represent the 3-year data and the asterisks 
the 1— year data. The bands represent the acceptance intervals as 
in previous figures. 



very similar results. The data lie outside the 1% accep- 
tance interval at scales Rg and i?io- The 3-year data 
show slightly higher values than the 1-year data at these 
scales. In particular, the upper tail probability for the 
1-year data was 0.56%, whereas for the 3-year data we 
obtain 0.38% at scale R 9 . 
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Fig. 5. — The left panel shows the cold Area in pixels, at threshold 4.0 versus the number of the scale. In the right panel the cold Area 
is represented versus the thresholds, while the scale is fixed at Rg. As in previous figures the asterisks represent the 1-year and the circles 
the 3— year data. The bands represent the acceptance intervals as in Figure 01 




Fig. 6. — Histogram of all biggest spots of the simulations at 
threshold 4.0 and scale Rg. The dashed vertical line represents the 
Spot in the 1— year data and the solid one represents the Spot in 
the 3— year data. 



3.3. Area 

We define the hot Area as the number of pixels above 
a given threshold v and the cold Area as the number 
of pixels below a given threshold —v. The threshold is 
given in units of the dispersion of the considered map. 

In C05 the total cold Area of the 1-year data was found 



to deviate from the Gaussian behaviour at scales Rs and 
i?9 and thresholds above 3.0 (see Figures 1 and 2 in C05). 

C05 found that the large cold Area of the Spot was 
responsible for this deviation. Such a big spot was very 
unlikely to be found under the Gaussian model at several 
thresholds (see Table 2 of C05). 

In the present paper we will define the Area as the 
maximum between hot and cold Area at a given threshold 
and scale. As for the Max estimator, we obtain in this 
way a more conservative estimator since the Spot will 
be compared to the biggest spot in each simulation no 
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Fig. 7. — Higher Criticism values of the 1-year WCM (aster- 
isks) and the 3-year WCM (circles). The acceptance intervals are 
plotted as in previous figures. 



matter if it is a cold or a hot spot. 

However the Area still deviates from the Gaussian be- 
haviour as can be seen in Figure The most significant 
deviation is again found at scale Rg and thresholds above 
3.0. 

Figure El shows the histogram of the biggest spot of 
each simulation compared to the 1-year and 3-year Area 
of the Spot at scale Rg and threshold 4.0. The Spot is 
more prominent in the 3-year data and only very few 
simulations show bigger spots. The upper tail proba- 
bilities obtained at scale Rg for 1-year and 3-year data 
are presented in table |21 As in the previous estimators, 
the 3-year data are in general slightly more significant. 
The new and more conservative definition of the Area 
estimator reduces the upper tail probability of the Spot 
although it is still widely below 1%. 

3.4. Volume 

From the previous subsections we know that the Spot 
is extremely cold and it has a large Area at thresholds 
above 3.0. The best estimator to characterise the Spot 
would be therefore the volume. Hence we define the 
volume referred to a particular threshold as the sum of 
the temperatures of the pixels conforming a spot at this 
threshold. In Table|3]we compare the probability of find- 
ing a spot with higher or equal Volume as the data, as- 
suming the Gaussian hypothesis. The values are very 
similar to those obtained for the Area estimator. Val- 
ues for the Volume are slightly more significant and they 
show less variations with the threshold. 

3.5. Higher Criticism 

The HC statistic proposed by Donoho & Jin (2004) was 
designed to detect deviations from Gaussianity that are 
caused by either a few extreme observations or a small 



Fig. 8. — Higher Criticism of the 3-year WCM at scale Rg. 
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Fig. 9. — Image projected as in Figure HI showing the 3-year 
WCM map (upper panel) and the Higher Criticism map (lower 
panel), both at scale Rg. 



proportion of moderately extreme observations. More- 
over, the statistic provides a direct method to locate 
these extreme observations by means of HC values cal- 
culated at every individual data point. 

For a set of n individual observations Xi from a cer- 
tain distribution (Xi normalized to zero mean and dis- 
persion one), HC is defined as follows. The ob- 
served values are first converted into p-values: pa\ — 
P{\N(0, 1)| > \Xi\}. After sorting the p-values in as- 
cending order < pp) < • • • < P( n ), we define the HC 



at each pixel with p- value pi, by: 

i/n - P(i) 



HC n . 



TABLE 4 

p-VALUES FOR DIFFERENT ESTIMATORS. 



We compute the values of the HC statistic of the 3-year 
WCM in real and in wavelet space. The obtained values 
of the HC statistic are presented in Figure These val- 
ues correspond to the maximum of the HC values found 
at the individual pixels. As in previous figures, circles 
denote the results obtained from the 3-year WCM, as- 
terisks those from the 1-year WCM and the bands repre- 
sent the acceptance intervals. As one can see in the Fig- 
ure, the data in wavelet space are not compatible with 
Gaussian predictions at scales R$ and Rg at the 99% c.l. 
This is in agreement with the result obtained by CJT for 
the 1-year WMAP data although there the HC values 
at scale i?g were just below the 99% c.l. The upper tail 
probabilities for the 1-year and 3-year maximum HC val- 
ues at scale Rg, are 0.56% and 0.36% respectively. The 
map of HC values at scale Rg is presented in Figure |5| It 
is clear that the pixels responsible for the detected devi- 
ation from Gaussianity are located at the position of the 
Spot. Convolution with the wavelet causes the observed 
ring structure in the HC map. Figure El shows a blowout 
image of the Spot as it appears at scale Rg in the wavelet 
map and in the HC map. 

4. SIGNIFICANCE 

In the previous section, the upper tail probabilities of 
each estimator at scale Rg were given. All the consid- 
ered estimators showed the lowest upper tail probability 
at scale Rg. However these are not rigorous measures 
of the significance of the Spot, since the number of per- 
formed tests is not taken into account. In this section we 
will recalculate the p- value of the deviation in the kur- 
tosis found by V04 and discuss the issue of a posteriori 
significances. 

When an anomaly is detected in a data set following a 
blind approach, usually several additional tests are per- 
formed afterwards to further characterize the anomaly. 
In most of these cases, the only reason these tests have 
been performed is the previous finding of the initial 
anomaly. If another anomaly would have been detected, 
other followup tests would have been performed. Hence 
these followup tests have not been performed blindly and 
should not be taken into account to calculate the signif- 
icance of the initial detection. 

This issue was already discussed in C06 and McEwen et 
al. (2005). Both papers recalculated the significance of 
the excess of kurtosis in the 1-year WCM found by V04. 
The excess of kurtosis was found performing a blind test, 
since no model was used and no previous findings con- 
ditioned the choice of the scales. Since 15 wavelet scales 
and two estimators (skewness and kurtosis) were consid- 
ered, a total sum of 30 tests were performed. Three of 
these tests detected a strong deviation from Gaussianity. 
Scales i?7, i?8 and Rg presented upper tail probabilities 
0.67%, 0.40% and 0.38% in the 1-year data. This fact 
was taken into account in C06, but it was not by McEwen 
et al. (2005). The latter searched through the simula- 
tions in order to find how many of them showed a higher 
or equal deviation than the maximum deviation of the 
data, ignoring that the data showed a high deviation at 
two adjacent scales. The p-value found in this way was 



Estimators 


p- value 


kurtosis 


0.86% 


skewness + kurtosis 


1.85% 


Max 


11.64% 


Area 3.0 


3.27% 


Area 4.0 


1.09% 


Higher Criticism 


3.48% 



4.97% whereas C06 obtained 1.91% taking into account 
that the data deviate at three consecutive scales. It is 
also interesting to note that, when both Galactic hemi- 
spheres were considered independently, C06 found a p- 
value of 0.69%, although this could be considered as a 
followup test. 

Some readers could find that the three-consecutive- 
scales criterion is an a posteriori choice since we look first 
at the data and given that they deviate at three consec- 
utive scales, we then calculate from the simulations how 
probable this is. Therefore we should consider a new test 
which eliminates this a posteriori choice. We fix a priori 
a significance level which is the 1% acceptance interval 
given in all figures, and count for each estimator (skew- 
ness and kurtosis) how many scales lie outside, no matter 
if they are consecutive or not. Then we search through 
the simulations how many show at least that many scales 
outside the 1% acceptance interval as the data. 

Applying this test to the 3-year WCM, we find that 
scales R$ and Rg lie outside the 1% acceptance interval 
and scale R-? lies on the border for the kurtosis estimator 
as can be seen in Figure |3 Searching through the sim- 
ulations how many deviate in three scales either in the 
skewness or in the kurtosis estimator, we find a p-value 
of 1.85%, which is still below the p- value obtained for the 
1-year data with the three-consecutive-scales criterion. 

As already discussed we should not include the fol- 
lowup tests in a rigorous significance analysis. However 
it is difficult to assess if some of these tests would have 
been performed or not without the first finding of V04. 
In fact, the area and maxima analyses are very intuitive 
and simple. If V04 had performed their blind analysis on 
those estimators instead of using skewness and kurtosis, 
then the significance would be different. We should dis- 
tinguish between those tests which are clearly followup 
tests, because the only reason they have been performed 
is the initial detection, and other tests which just have 
been performed after the initial detection, but could have 
been performed before. 

Hence we apply our new robustness test to kurtosis, 
Max, Area at thresholds 3.0 and 4.0 and Higher Crit- 
icism separately. Note that whereas the first two esti- 
mators are two-sided, the Area and Higher Criticism are 
one sided estimators. The p-values obtained in this way 
are listed in Tabled The kurtosis and Area at thresh- 
old 4.0 show p- values around 1%, Higher Criticism and 
Area at threshold 3.0 around 3%. On the contrary the 
Max estimator does not show a significant deviation from 
Gaussianity according to this robustness test. 

The most conservative and reliable value is the 1.85% 
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Fig. 10. — Frequency dependence of the temperature at the 
center of the Spot at scale Rg. Again the asterisks represent the 1— 
year data and the circles the 3— year data. The horizontal line shows 
the value of the 3-year WCM. The data at the same frequency have 
been slightly offset in abscissa for readability. 
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Fig. 11. — Maximum absolute observation for the Q, V and W 
bands, compared to the 3— year WCM values. 



Fig. 12. — Frequency dependence of the Area of the Spot at 
scale iig and several thresholds. Asterisks represent the 1-year 
data and the circles the 3— year data. The 3-year WCM values are 
represented by horizontal lines. 
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Fig. 13. — Higher Criticism values for the Q, V and W bands, 
compared to the 3— year WCM values. 



figure since it is not suspicious of being obtained through 
a posteriori analyses. Nevertheless it is still noticeable 
that the followup tests performed in C05, C06, CJT and 
in the present paper, confirm the initial finding with a 
very similar significance. Even if strictly speaking these 
should not be taken into account for establishing the sig- 
nificance of the Spot, they confirm the robustness of the 
detection. 



5. FREQUENCY DEPENDENCE 

In this section we will analyse the frequency depen- 
dence of the previously analysed estimators. A flat fre- 
quency dependence is characteristic of CMB, whereas 
other emissions such as Galactic foregrounds show a 
strong frequency dependence. Figure IT~TI shows that the 
kurtosis has almost identical values at the three fore- 
ground cleaned channels, namely Q,V and W. Same be- 
haviour was observed in the 1-year data (see Figure 7 in 
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Fig. 14. — Kurtosis values for the Q, V and W bands, compared 
to the 3-year WCM values. 



issues raised by Spergel et al. (2006), since it is very dif- 
ficult to get completely rid of it. Most analyses perform 
many tests and it is not easy to assess how many of them 
are followup tests and which is the probability of finding 
an anomaly by chance. As discussed in section 4 a very 
careful analysis shows that the Spot remains statistically 
significant at least at the 98% confidence level, without 
using any a posteriori statistics. 

In addition C06 proved that the Spot remained highly 
significant no matter which foreground reduction tech- 
nique was used. These results are confirmed in the 
present paper. The new foreground reduction used in 
the 3-year data enhances slightly the significance of our 
detection. Moreover the multi-frequency analysis of the 
previous section shows an even flatter frequency depen- 
dence of the Spot. 

As already discussed in previous sections the noise does 
not affect significantly our wavelet analysis. In fact the 
coadded 3-year results are very similar to those obtained 
with the 1-year data of the new data release. No signif- 
icant cold spot is observed based on the analysis of the 
three difference maps (year 1 - year 2, year 2 - year 3, and 
year 1 - year 3). Moreover Figure ITTTI shows that even the 
particularly 1/f contaminated W4 Difference Assembly 
shows almost the same result as all the other Difference 
Assemblies. 



C06). Strong frequency dependent foreground emissions 
are unlikely to produce the detected excess of kurtosis. 

The frequency dependence of the temperature at the 
center of the Spot, i.e. at the pixel where the temperature 
of the Spot is minimum in the WCM map, is presented in 
Figure EH The error bars of the 1-year data have been 
estimated performing 1000 noise simulations as explained 
in section 5.1 of C06. As the noise variance is « 3 times 
lower in the 3-year data, we estimate the new error bars 
simply by dividing the old ones by No frequency de- 
pendence is found for the new data set in agreement with 
the results for the 1-year data. Max, Area and HC val- 
ues at different frequencies (see Figure ^2 Figure Upland 
Figure IT3"|) show a very low relative variation compared 
to the 3-year WCM. 

All these results confirm the analysis performed in sec- 
tion 5 of C06 where the data were found to fit a flat 
CMB spectrum. The present analysis confirms the dis- 
agreement between the conclusions of C06 and those of 
the work of Liu & Zhang (2005) where Galactic fore- 
grounds were considered to be the most likely source for 
non-Gaussian features found with spherical wavelets. 

6. DISCUSSION 

Spergel et al. (2006) enumerate several reasons to 
be cautious about the different anomalies found in the 
WMAP data: Galactic foregrounds or noise could be 
generating the non-Gaussianity, and in addition most of 
the claimed detections are based on a posteriori statis- 
tics. Also spatial variations of the noise variance and 1/f 
noise could affect some of the perfomed analyses. They 
suggest several tests to be done using difference maps 
(year 1 - year 2, year 2 - year 3, etc.) and multi-frequency 
data. 

We have tried to address all those points for the Spot. 
The a posteriori analysis is one of the most important 



7. CONCLUSIONS 

In this paper we repeat the analyses that detected the 
non-Gaussian cold spot called the Spot at (b = —57°, I — 
209°) in wavelet space in the 1-year of WMAP data, us- 
ing the recently released 3-year WMAP data. The pre- 
vious works V04, C05, CJT and C06 found the Spot to 
deviate significantly from the Gaussian behaviour. The 
Spot was detected using several estimators, namely kur- 
tosis, Area, Max and HC. This work confirms the de- 
tection applying all these estimators to the recently pub- 
lished 3-year WMAP data. At scale i?g, the upper tail 
probabilities of all these estimators when applied to the 
3 year WMAP data are smaller than the corresponding 
ones for the first year WMAP data. This is mostly due 
to the improved foreground reduction of the data. We 
calculate the probability of finding such a deviation from 
Gaussianity considering only skewness and kurtosis since 
these were initially used by V04 following a blind ap- 
proach. Therefore excluding followup tests which could 
be considered as a posteriori analyses we obtain a p- value 
of 1.85%. Moreover, the Spot appears to be almost fre- 
quency independent. This result reinforces the previous 
foreground analyses performed by C06. It is very unlikely 
that foregrounds are responsible for the non-Gaussian be- 
haviour of the Spot. Comparing the WMAP single year 
sky maps, we conclude that the noise has a very low con- 
tribution to our wavelet analysis as already claimed in 
V04, C05. Future works will be aimed at finding the ori- 
gin of the Spot. As discussed in the introduction several 
possibilities have been considered, based on Rees-Sciama 
effects (Rees & Sciama 1968, Martinez-Gonzalez & Sanz 
1990, Martinez-Gonzalez et al. 1990) and inhomogenous 
or anisotropic universes. We are presently working on 
studying another: topological defects (Turok & Spergel 
1990, Durrer et al. 1999) as textures could produce cold 
spots. New and more detailed analyses arc required in 
order to answer that question. 
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